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ABSTRACT 


It was suggested by Munk and Forbes (1989) that climate induced changes in ocean 
temperature may be monitored by measurements of cross-basin acoustic travel time 
variability. The feasibility of such a monitoring system depends on the spatial and 
temporal variability of the cross-basin acoustic paths in the presence of ocean variability 
of many different scales. For this thesis the variations in arrival position, azimuthal ar¬ 
rival angle, ray trajectory and the corresponding changes in travel times along the 
three-dimensional multipaths due to meso- and gyre scale ocean temperature fluctu¬ 
ations were analyzed. Emphasis was placed on the acoustic paths from Heard Island in 
the Indian Ocean, the proposed location for the sound source, to the west coast of the 
United States. 

It was found that the launch angles of reliable, unimpeded acoustic paths to the west 
coast of the United States lie within a 3° x 2° angular sector. The variability of this rav 
envelope is discussed. 

An optimal receiver site location was found to exist in the vicinity of Monterey Bav. 
California. The possibility of a proposed listening site location near Coo* Bay, Oregon, 
was also examined. However, the ray paths to Coos Bay interact with the bottom fre¬ 
quently, thus rendering them less reliable. 

All the ray traces for this study were carried out using the recently upgraded 
Hamiltonian raytracing code HARPO, interfaced with output from the Semtner-Chervm 
eddy resolving global ocean general circulation model. This interface allows for a real¬ 
istic simulation of the effects of meso- and gyre scale processes on the variability of the 
various cross-basin paths. The following alternative interfacing methods were investi¬ 
gated: Hardy's multiquadric method and an empirical orthogonal functions method 
developed by Newhall, el al., 1989. 
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I. INTRODUCTION 


A. THE GREENHOUSE EFFECT 

Investigation of potential global warming trends, caused by increasing 
atmospheric content of gases like carbon-dioxide, methane, 
chlorofluorocarbons (CFC's) and of other trace constituents, is one of the 
most challenging research topics for this and future decades. Basically, 
these so called “greenhouse gases” trap heat in the layer between the 
earth's surface and the level from which radiation escapes back to outer 
space by absorption of the long wave back radiation from the earth's 
surface. Increased levels of these gases in the atmosphere will thus raise 
the earth's temperature, i.e., the “greenhouse effect”. Almost everybody 
agrees on the fact that greenhouse gas concentrations in the atmosphere 
have increased significantly. Since 1765, ievels of atmospheric carbon 
dioxide have increased from 280 parts per million by volume to more than 
350 parts per million. Methane concentration has more than doubled, 
from 800 to 1700 parts per billion, and nitrous oxide has increased by 
about ten percent, from 285 to 310 parts per billion. During the past 30 
years, concentrations of CFC's have risen from essentially zero to about 
one part per billion (Jones and Wigley, 1990). However, much contro¬ 
versy exists over related questions concerning the amount of global 
warming these increases will cause, the time span over which it will take 
place and whether trends can be discerned and possibly extrapolated to 
the future. Actual observations of atmospheric surface warming suffer 
from the fact that many of the land stations are contaminated by signif¬ 
icant changes in the micro-climate, especially in urban areas, while sea 





surface temperature measurements are biased by a conversion from 
bucket to engine cooling-water intake measurements. Additionally, it is 
difficult to measure greenhouse warming in the atmosphere, because it 
has similar spatial patterns as the background air temperature variability. 
In the ocean the natural temperature variability is significantly different 
from the expected greenhouse signal. Accordingly, the oceans may be a 
better environment for early detection of greenhouse warming (Barnett 
and Schlesinger, 1987). 

B. THE HEARD-ISLAND FEASIBILITY EXPERIMENT 

It was suggested by Munk and Forbes (1989) that climate induced 
changes in ocean temperature can be monitored by measurements of 
cross-basin acoustic travel time variability. The possibility of detection 
of acoustic transmissions over cross-basin ranges was demonstrated bv a 
1960 experiment, in which detonations at the depth of the sound channel 
axis off Perth, Australia, were recorded on hydrophones located near the 
depth of the sound channel axis in the vicinity of Bermuda (Shockley. 
et al., 1982). In addition, a bottom -mounted hydrophone in the same 
location also detected significant signal. This result showed that, apart 
from propagation along the sound channel axis, off-axis cross-basin 
propagation could also survive and that experimental data could contain 
information about the temperature field in layers other than only the 
axial one. 

Munk and Forbes (1989) estimated that greenhouse effects warm the 
sound channel axis at a rate of 0.005 °C per year, which translates into a 
decrease of acoustic travel time on the order of 0.1 s per year along the 
cross-basin axial paths. Since meso- and gyre scale processes generate 
about 1 s rms fluctuations in the cross-basin acoustic travel times. 




measurements must l e carried out over at least a decade. A field 
experiment led by the Applied Physics Laboratory in Seattle, Washington, 
is planned for January 1991 to investigate the feasibility of such a decadal 
monitoring system. The plan provides for a hydroacoustic low frequency 
source array to be lowered to about i 14 m depth in 1000 m of water 
southeast of Heard Island (53.2°S / 73.1°E), off Australia in the Indian 
Ocean. The sources are to be powered from the attending vessel, R V 
Cory Chouest. Initially 12 days of transmissions are planned, cycling 
through four types of signals: continuous-wave at 5 7 Hz, a five-tone-like 
signal called pentaline. a 255 digit maximal-linear-shift register phase- 
coded sequence, commonly used for tomographic experiments (Spindel, 
1985) and longer m-sequences signals. Transmissions will commence on 
January 24, 1991. The source site near Heard Island was chosen, because 
it has access to all the major ocean basins. It allows for unimpeded re¬ 
fracted ray paths to Bermuda, the west coast of the United States and 
possibly a path through the Tasman Sea reaching Coos Bay, Oregon, as 
can be seen in Figure 1. 

Investigators from the U.S. Naval Postgraduate School, the Monterey 
Bay Aquarium Research Institute and the Massachusetts Institute of 
Technology will jointly deploy a listening array in the vicinity of 
Monterey Bay on the west coast of the United States. 

The feasibility experiment addresses several major uncertainties 
underlying global acoustic transmissions: 

• What is the required source strength for a permanent installation and 
what type of acoustic source is optimal for this purpose ? 

• How is the variability of cross-basin acoustic ray paths affected by 
spatial ocean variability ? 
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Figure 1. Cross-basin acoustic rays from Heard Island 
Source: Munk and Forbes 19S9. 


• Will the arrival pattern of the acoustic signals have sharp and stable 
features that allow for determination of the variable travel times with 
at least 10 ins accuracy ? 

• What are the optimal locations for the listening stations and what is 
the best type of receiver array ? 

If the 1991 experiment is successful in demonstrating the feasibility of 
global acoustic transmissions by resolving these uncertainties, a perma¬ 
nent climate monitoring system could be installed. This permanent sys¬ 
tem will provide early warning of upper ocean warming in the southern 
oceans and of interior warming in other oceans. 
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In addition to this possibility of monitoring global ocean warming, a 
permanent cross-basin range acoustic transmission system could enhance 
our capability in predicting the ocean's meso- and gyre scale features. 
Travel time variabiliu -dong paths going through different ocean basins 
would contain information on meso- and gyre ale processes in those 
basins, that can be used for the calibration of eddy resolving global ocean 
circulation models. The statistics of the measured travel time variability 
can be used to check the statistics of the global, general ocean circulation 
model output from time to time for consistency. 

C. THESIS OBJECTIVES 

This thesis is part of ongoing computer simulation studies of cross- 
basin acoustic transmissions in support of the analysis and modeling of 
the field data that will be obtained during the Heard Island feasibility 
experiment in January 1991. The effects of ocean fronts and eddies on 
the acoustic ray paths crossing the Indian and Pacific Oceans to the west 
coast of the United States from a source near Heard Island are investi¬ 
gated. Specifically, the primary objectives of this thesis research are: 

* to determine the locations along the west coast of the United States, 
which will be insonified by the Heard Island sound source; 

• to examine and, where possible, quantify spatial and temporal 
variations in arrival position, azimuthal arrival angle, ray trajectory 
and the corresponding fluctuations in travel time along the three 
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dimensional acoustic multipaths, due to meso- and gyre scale ocean 
temperature fluctuations; 

• to investigate whether an optimal listening site location exists in the 
vicinity of Monterey Bay, California; 

• to investigate the potential shadowing effects on the arriving acoustic 
rays of the shallow “Erben” and “Fieberling” seamounts, located off 
the California coast; 

• to examine the possibility of unimpeded reliable ray paths from the 
proposed Heard Island source location through the Tasman Sea to 
Coos Bay, Oregon. 


Using numerical models, this thesis attempts to answer some of the cru¬ 
cial questions concerning the usefulness of a global acoustic monitoring 
system and whether it can be realized. The results presented here will aid 
in the interpretation and analysis of the experimental data to be obtained 
in the feasibility field experiment of January 1991. Results will have 
definite implications for the design of a permanent global acoustic trans¬ 
mission experiment, in terms of listening site location, expectable arrival 
structure, etc.. 

Naval operations can benefit from this research as well. Results will 
provide a better understanding of the effects of mesoscale processes on 
cross-basin low frequency sound propagation in support of naval oper¬ 
ations. Western navies have an increasing interest in low frequency long 
range acoustic prediction research in support of antisubmarine warfare 
efforts. Low frequency sonar operations, both active and passive, will 
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benefit from greater knowledge of the characteristics of low frequency 
three-dimensional vvavefields that propagate through oceanic fronts and 
eddies, while integrated acoustic and physical oceanic modeling capabili¬ 
ties in this field are of prime importance to these operations. 

An accurate long range acoustic monitoring and prediction system can 
only be established through the integration of physical ocean modeling 
efforts with three-dimensional acoustic modeling. In this thesis the results 
of such an integration are presented. Prior work in this subject area has 
only involved the coupling of acoustic models to regional, small domain 
ocean models (Chiu and Ehret, 1990). The following alternative interpo¬ 
lation methods for the interfacing procedure to link the output data sets 
from the Semtner-Chervin global ocean general circulation model to the 

three-dimensional acoustic raytracing model are examined: 

• Hardy's multiquadric interpolation scheme (Hardy, 1990) modified 
for three-dimensional earth-centered spherical coordinates; 

• the Empirical Orthogonal Functions (EOF) method, developed by 
Newhall, et al., 1989. 

The comparison between the two methods will aid the continuous effort 
to examine whether existing procedures and methods in the ongoing 
acoustic computer simulation studies can be improved in terms of accu¬ 
racy and efficiency. 
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D. THESIS APPROACH AND OUTLINE 


A three-dimensional Hamiltonian raytracing code, developed ori¬ 
ginally by Jones, R.M. et al. (1986) for a CDC-computer, HARPO, was 
selected to perform the raytracing in this thesis. Global output data sets, 
with mesoscale resolution in both space and time, generated by the 
Semtner-Chervin eddy resolving global ocean general circulation model 
(1988) were interfaced with the HARPO raytracing code. Simulations 
were carried out on a UNIX Sun workstation, with a 32 Mbyte random 
access memory and a 640 Mbyte hard disk memory. In Chapter II the 
theoretical background of both the ocean- and the acoustic models will 
be presented. It will contain a discussion on Hamiltonian raytracing and 
on the global ocean eddy resolving general circulation model. Two al¬ 
ternative methods for the interpolation of the input sound speed fields to 
the acoustic model from gridded temperature and salinity output data of 
the Semtner-Chervin ocean model will be discussed in Chapter III. The 
first method uses spiine fits to vertical Empiiical Orthogonal Functions 
(EOF's) and EOF coefficients and has been shown to perform well for 
small regional data sets (Newhall, et al ., 1989). Secondly the performance 
of Hardy's multiquadric technique, modified for three-dimensional 
earth-centered spherical coordinates, was tested against observed 
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conductivity, temperature and depth (CTD) data, acquired during the 
operational oceanography cruise on board RjV Point Sur from 5 to 9 
May 1990. Based on these tests, the multiquadric technique was rejected 
and the interfacing procedure, developed by Newhall, et al., (1989), was 
used for the raytracing simulations. 

Raytracing was performed using the HARPO code for a 900-day mean 
field ocean model output data set, as well as 11 consecutive instantaneous 
fields, extracted from the model output at 30-day simulation time inter¬ 
vals. Rays were traced from the Heard Island source location, at 1° in¬ 
terval launch azimuth angles and a 2° downward elevation launch angle. 
All the rays that made it through the Indian Ocean, passing south of New 
Zealand, and the Pacific Ocean to the west coast of the United States 
were used in the ray variability analysis. Additionally, the possibility of 
unimpeded rays through the Tasman Sea, north of New Zealand, was 
investigated. The results of these analyses are presented in Chapter IV. 
Variations in arrival position, azimuthal arrival angle, and ray trajectory, 
as well as the corresponding fluctuations in travel time along the three- 
dimensional multipaths were estimated. These estimates along with esti¬ 
mates of the locations along the west coast of the United States which 
will be insonified by the Heard Island sound source and of an optimal 
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listening site location in the vicinity of Monterey Bay are also presented 
in this chapter. In Chapter V a summary of this study and its conclusions 
is presented, along with some recommendations for improving future 
modeling work of this kind. 
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II. THEORETICAL BACKGROUND 


A. HAMILTONIAN ACOUSTIC RAYTRACING 

Many practical problems in ocean acoustics can be solved by applying 
geometrical acoustics, or ray theory. This technique allows one to simu¬ 
late the propagation of sound waves through a medium whose refractive 
index structure varies in a complex manner. The calculations of some 
acoustic quantities cannot be easily made in any other way, for example, 
computing multipath pulse travel times or showing which parts of the 
ocean affect each pulse arrival. However, many raytracing models, which 
use layers or cells where each raypath segment can be computed in closed 
form, do not take full advantage of the power of geometrical acoustics. 
Discontinuous sound velocity gradients at layer interfaces may cause false 
caustics to occur and cause discontinuous behavior of ray properties as 
launch angles vary. Often it is also difficult to extend these models to 
three-dimensional oceans, to account for currents and to compute re¬ 
flections off complicated bottom topography (Jones, R.M., et al., 1986). 

The acoustic model, which is used for this study, is the Hamiltonian 
Raytracing Program for the Ocean, HARPO. It was designed to over¬ 
come the above mentioned limitations. The user defines an ocean model 





by writing closed-form expressions for its sound speed and current dis¬ 
tribution in three dimensions, and by defining the bathymetry as a func¬ 
tion of geographical position. Because it uses continuous models, 
HARPO avoids the false caustics and discontinuous raypath properties 
encountered in conventional raytracing methods. 

In Hamiltonian raytracing, each raypath is calculated by numerically 
integrating Hamilton's equations with a different set of initial conditions. 
These equations are written in terms of a Hamiltonian function, 
H{p x ,p 2 , q x , q 2 ,....,q n ) , which is the system's total energy (kinetic plus 

potential), expressed as a function of generalized coordinates </,, q 2> . q n , 

and generalized momenta p x ,p 2 , Hamilton's equations have the 

following form: 


d Pi _ dH 
dt 8q i 


i = 1,2,—,/z. 


( 1 ) 


dc h _ dH 

dt dp i 


i = 1,2,—,/z. 


( 2 ) 


Lighthill (1987) derived Hamilton's equations for acoustic waves in four 
dimensions (three spatial and one temporal) for cartesian coordinates, 
using the “correspondence principle" from quantum mechanics. The 
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three cartesian coordinates x t for a “pocket” of wave energy travelling 
along its ray, correspond to the generalized coordinates q t in the dynam¬ 
ical system, the wavenumber /c, corresponds to the momentum p, and the 
radian wave frequency co corresponds to the energy H. Lighthill further 
showed, that co is constant along rays, which is equivalent to the fact that 
for any solution of the Hamiltonian equations of motion the total energy 
H remains constant. This results in the following system of “ray- 
equations”: 


dkj 

~li 


da) 

OX; 


i= 1,2,3. 


(3) 


dxj 

~dt 


dco 
dki 


i= 1,2,3. 


(4) 


Converting these equations to earth-centered spherical polar coordinates, 
where r is the radius, 4> the longitude and 6 the co-latitude, gives the 
following system of six first order equations that are used in the HARPO 
raytracing code (Jones, R.M., et al., 1986): 


dr _ dH 
di dk r 


( 5 ) 
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( 6 ) 


dd 1 oH 

dr r dkQ 

d<t> __1 8H 

dr r sin 6 dk^ 


(7) 


dk r 

dr 


( 8 ) 


dk 6 

dr r y 


cH 

86 



+ k^r cos 6 



(9) 


dk 4> 

dr 


—!— ( 
r sin 6 \ 


8H , ■ u dr , a d6 

~r~ - kj, sin 6 —- k^r cos 6 -j- 

c<p 0 dr V dr 


( 10 ) 


Here r is the independent variable, sometimes proportional to time, 
whose physical meaning depends on how the Hamiltonian is defined. 
For acoustic waves in the ocean, the Hamiltonian, which is constant along 
a raypath, is defined as: 

H(R,k) = [v-k (r) • V - c 2 (r) k 2 (r) = 0 , (11) 

where is the ocean current, c(jij is the sound speed Field, and co is 

the angular wave frequency. V and c may also depend on time. Thus, the 


14 





effects of a three dimensional vector current field are automatically in¬ 
cluded in the definition of the Hamiltonian. 

A few final remarks must be made here about the HARPO code: 

• Its computations lie entirely within the scope of geometrical acous¬ 
tics. It does not apply corrections for diffraction or partial reflections. 
The ocean model must be deterministic, not random. 

• It does not compute acoustic intensity directly. Also, it is not set up 
to find eigenrays, that connect source and receiver. However, it is 
relatively straightforward to write supplementary code to carry out 
these computations. 

• It cannot handle refraction at discontinuities of refractive index or 
its gradients. 

• It does not make checks to see if the ocean model satisfies conserva¬ 
tion laws and boundary conditions, or that current and sound speed 
models are geostrophically consistent. 

• It integrates the system of Hamiltonian ray equations using the 
Adams-Moulton predictor-corrector method with a Runga-Kutta 
starter. Initial conditions are specified in terms of ray launch angles 
(wavenumber components). 

• Since it works in spherical polar coordinates, earth curvature effects 
are accounted for automatically. However, effects of the earth's 
ellipticity that cause geodesic veering of the acoustic rays are ignored. 

B. THE GLOBAL OCEAN DATA SETS 

The continuous ocean model required in the HARPO program is 
constructed from Semtner-Chervin eddy resolving global ocean general 
circulation model (SC-model) output data sets. The output of this ocean 
model was chosen, because it constitutes the only global data set available 
today, that has mesoscale resolution spatially as well as temporally and 
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because it has been extensively tested with satisfactory results (Chiu, 
et a /., 1990). The SC-model evolved from the Bryan and Cox model 
(Bryan and Cox, 1967). Through efforts by several different investigators, 
this original code was revised to allow for arbitrary bottom topography 
and coastlines. Also, it was updated for vector-processing and for 
supercomputers employing multiple parallel processors (Semtner and 
Chervin, 1988). 

The SC-model is a primitive equation, free-thermocline global ocean 
circulation model. Its domain extends from 75°S to 65°N and through 
360° of longitude, with cyclic continuity at 20°E. It has 1/2° x 1.2° grid 
spacing in the horizontal and 20 unevenly spaced levels in the vertical. 
To obtain ocean depths for the model, a data set of global bathymetry 
available at the National Center for Atmospheric Research for 1/12° lat¬ 
itude and longitude points was spatially averaged over 1/2° x 1/2° boxes 
centered at the model gridpoints. The bathymetry data was smoothed 
twice, to eliminate topographically induced numerical instability. This 
procedure resulted in an average ocean model depth of 4400 m and a 
maximum depth of 5200 m. In addition, some simplifications in geom¬ 
etry were made to create a global 1/2° data set of continental outlines. 
The Arctic was walled off at 65°N latitude. The Red Sea, the Gulf of 
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Aden, the Persian Gulf, the Gulf of Oman, the Andaman Sea and the 


Great Bight of Australia were all filled in. Some major islands were 
connected to the continental land masses: Madagascar was attached to 
Africa, the main Indonesian islands were incorporated into a larger, 
simplified southeast Asian peninsula and the islands of New Guinea and 
Tasmania were connected to Australia. Other, smaller, islands were re¬ 
placed with shoals of 100 m depth. This was done, to allow their major 
topographic features to remain in place, without the need for elaborate, 
rigorous numerical treatment of these islands. However, the replacement 
of certain islands by shoals might prove to be troublesome from an 
acoustic point of view. Depending on the ray trajectories, acoustic energy 
that would otherwise be impeded by these islands, could now possibly 
leak out through the introduced shoals. Since the ray paths expected in 
the set up of this analysis will traverse depths that are a lot deeper than 
100 m, it is unlikely that this problem will affect the results. 

Three “real” islands remained, namely. New Zealand, Australia-New 
Guinea and Antarctica. The simplified model geometry is presented in 
Figure 2. 
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Figure 2. The Semtncr-Chervin global ocean model geometry 

Source: Semtner and Chervin, 1988. 

In Figure 3 contours of the resulting bathymetry from 500 to 3000 m 
for the area of interest in this thesis research are presented. The contour 
interval is 500 m, and not all of the geometry simplifications are shown. 

It should be noted that two shallow areas exist, that might impede 
acoustic rays. The first area south of New Zealand is called the Campbell 
Plateau. The second area is the Tasman Sea, between New Zealand and 
Australia. Additionally, there exist two shallow seamounts off the 
California coast, that were omitted from the modeled bathymetry. These 
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Figure 3. Bathymetry from the Semtner-Chervin ocean model: Contours from 
500 to 3000 m, contour interval 500 m. 

are: “Fieberling” in position 32.4°N / 23'\7°E, with a depth of 431 m 
and “Erben” in position 32.9°N / 227.5°E, with a depth of 921 in. In 
Chapter IV it will be examined whether the acoustic ray trajectories could 
possibly be affected by these seamounts. 

T he 20 vertical levels of the model are unevenly spaced, in a way to 
give enhanced resolution in the oceanic thermocline, with 12 levels in the 
upper 1400 m of the water column. Since acoustic ray paths generally 
will be confined to depths near the sound channel axis, that normally fall 
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within the upper 1400 m of the water column as well, the vertical structure 
of the ocean model will also be adequate from an acoustic point of view. 
The vertical “boxes” are stacked downward, with the local bathymetry 
determining the number of boxes. The vertical velocity is defined at the 
interface between boxes. T, S, u and v are all defined at the same vertical 
level. T and S are defined at the center of the boxes, while u and v are 
defined at the midpoints of the vertical edges of the boxes (Long, 1990). 
This vertical arrangement of grid points is illustrated in Figure 4. 

Forcing for the model consists of annual mean wind stress from the 
Hellerman and Rosenstein (19S3) global data sets, interpolated to the 
1,2° surface velocity gridpoints. A three-dimensional Field of temperature 
(T) and salinity (S) interpolated from the gridded data base of Levitus 
(1982) is used to force the model at the surface on a monthly time scale 
and at sub-thermocline depths on a three year time scale. The model was 
integrated for a simulated time of 22.5 years. Output data were archived 
at a simulation time interval of three days. With the application of the 
described forcing method, realistic meso- and gyre scale ocean variabil¬ 
ities resembling many of the observed features in the global ocean resulted 
from the simulation in the final model years. The mean field output data 
set, used in this thesis corresponds to a 900-day time average of the data 
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simulated between model years 20 and 22.5. In addition, 11 “monthly” 
fields, consisting of the instantaneous model output data at a simulation 
time interval of 30 days, were used to analyze meso- and gyre scale effects 
on the spatial and temporal acoustic ray variability. In the next chapter 
the interfacing procedure, which links these data sets to the acoustic 
Hamiltonian raytracing model, will be discussed. 
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III. OCEAN-ACOUSTICS INTERFACE 
A. INTRODUCTION 

The Semtner-Chervin eddy resolving global ocean general circulation 
model generated approximately 200 Mbyte of data every three simulation 
days. The data are archived at the National Center for Atmospheric 
Research. Because the acoustic raytracing simulations were to be carried 
out on a SUN-4 workstation at the U.S. Naval Postgraduate School, 
which has only 32 Mbyte of random access memory and a 640 Mbyte 
hard disk, the task to interface HARPO to the global ocean data set is 
not a trivial one. To get around this storage problem, two steps are taken 
(Chiu, et al., 1990): 

• Temperature and salinity data in individual ocean basins are first 
subsampled from the global data set, using a FORTRAN code de¬ 
veloped by Altman, 1990. These subsampled data values are then 
converted to sound speeds using the Mackenzie empirical formula 
(1981) which has a standard error of 0.070 m/s. Only the subsampled 
and thus converted sound speed data are transferred to the local 
SUN-4 workstation. 

• The subsampled sound speed data set is still considered large in a 
workstation environment. Processing all the data at the same time 
would cause the acoustic computations to come to a virtual halt. 
To make the simulations tractable, a rectangular window of 10° lati¬ 
tude by 20° longitude is used to move in a discrete fashion with the 
ray being traced. Only data inside a sub-region with boundaries de¬ 
fined by this window are extracted from the subsampled sound speed 
data set. 


23 




The interfacing procedure furthermore involves the selection of an effi¬ 
cient and relatively optimal interpolator for the generation of continuous 
sound speed fields, as required by the HARPO code. In this selection 
process, two different interpolation techniques were examined: 1) Hardy's 
multiquadric method modified for three-dimensional earth centered 
spherical, poiar coordinates and 2) an empirical orthogonal functions 
(EOF) method, developed by Newhall, et al., (1989). Both methods will 
be discussed here. 

B. HARDY'S MULTIQUADRIC METHOD 

The multiquadric (MQ) method was first discovered in 1968 (Hardy, 
1990). It is characterized in two dimensions by the following idea. For 
each position (x,-,y ( ) simply choose a basis function i?,(x,y) and then 
determine coefficients D t so that G {x ,y) = Y.d i B i {x ,y) interpolates the 

i 

data (Franke, 1979). The basis function of choice in Hardy's MQ scheme 
is the upper hyperboloid By = ^/(x; — x ; ) 2 -P {y t — y,) 2 + S 2 , where <5 is a 
parameter which determines the semi axis of the hyperbolas. Hardy's 
method can thus be described by the following formula: 


G(x,, y,) = dj y] (x; - xjf + 0/ - yj) 2 + d 2 i — 1,2(12) 

J= 1 
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Or in matrix notation: 


Gi = Bijdj (13) 

The coefficients are found by solving a set of linear equations in terms 
of the basis functions: 


<1, = B„ ' G, (14) 

and the evaluation or prediction (interpolation) at any point p is: 

Gp-Wi'G, ( 15 ) 

It must be noted, that for S = 0 the basis function in the interpolation 
scheme is reduced to a simple measure of distance away from grid points, 
which diverges linearly with increasing distance. However, non-zero val¬ 
ues for <5 must be used (even if they are infinitesimally small) to make 
each term in the scheme infinitely differentiable. For application in 
HARPO this is a necessary property of the interpolation method to be 
used. Theoretical investigations have proved that the MQ interpolation 
is always solvable for distinct data. Franke (1979) performed a critical 
comparison of about 30 interpolation methods with respect to their 
capability of reproducing several mathematical surfaces from sparse. 
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scattered data. He concluded after extensive testing, that Hardy's 
multiquadric method could perform very impressively. 

For this investigation the MQ interpolation scheme was modified for 
application to three-dimensional gridded sound speed data in spherical 
coordinates. Few people have applied MQ interpolation to more than 
two dimensions. However, theoretically the method can be applied to any 
number of dimensions without altering the scheme's basic properties. 
The basis function in the MQ scheme was transformed in terms of the 
spherical, polar coordinates r, 0 and 0: 
By = yj(r — r,) 2 + (dj — dj) 2 + (0, — 0 y ) 2 + S 2 . Equation (12) can now be 
rewritten as follows: 


n I - 

G ('i , 9/ , = £ dj 0r ; - rf + (0, - d/ + «>,- - 0/ + (16) 

for i = l,2,...,n. The coordinates r, 0 and d are defined as earth centered 
spherical, polar coordinates, where r is the radius, 0 is the longitude and 
0 is the co-latitude. This modification also allows for easy application 
of a correction for the earth's ellipticity. One could simply define the 
local radius, r, as a function of latitude in accordance with a given 
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ellipsoid of choice, using standard geodetic formulas. For cross-basin 
range acoustic transmissions, this correction is probably not negligible. 

To improve the accuracy of the interpolation and to obtain consistent 
results, the independent coordinates were scaled by dividing the respective 
terms by appropriate scale parameters. The term involving the radius was 
scaled by the square of the average ocean depth in kilometers, the longi¬ 
tude term was scaled by the square of the east/west width of the area in 
radians and the latitude term was scaled by the square of the north south 
width of the area in radians. This scaling effectively normalized each of 
the terms to values between zero and one. 

A problem that couia be encountered involves so called “track data"; 
i.e., data points which are closely spaced along one coordinate direction 
and widely spaced along the two orthogonal directions. Among others, 
Kansa (1990) found that “track data” give MQ interpolants with very 
large errors. In our setup we also had to consider “track data”, since the 
unsealed r coordinates have much larger values than the unsealed hori¬ 
zontal coordinates 0 and 8. However, from existing knowledge of ocean 
physics it is clear that aside from the uneven spacing of coordinates, we 
also have to allow for the horizontally stratified nature of the ocean. In 
effect this causes distinct sound speed values that reside at one and the 
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same vertical level to be much higher correlated than distinct sound speed 
values that reside at different vertical levels. To solve for this problem, 
the normalized integration formula was weighted in such a way, that 
components in the horizontal <j) and 6 directions seem much smaller than 
the distance component in the vertical r direction. This assigned a higher 
“importance” to the sound speed values at grid points that reside at the 
same vertical level as the position, for which we want to interpolate. At 
the same time this weighting solved the “track data” problem, as de¬ 
scribed above. Appropriate values for weighting the terms in the 
interpolation formula were determined experimentally. 

The modified MQ interpolation scheme was tested using data ac¬ 
quired during the operational oceanography cruise, from 5 to 8 May 1990. 
During this cruise, CTD casts were carried out from the R!V Point Sur 
at regularly spaced stations in the area between 36.0°-36.7°N and 
121.0°-123.0°W, off the California coast. From the CTD data at these 
stations sound speed values at 20 depth levels were computed from tem¬ 
perature, salinity and depth data, using the Chen and Millero (1977) 
sound speed formula. This data constituted the “true” observed data. 
The MQ method was tested by using one part of this data and recovering 
the other part by interpolation. Interpolation results were examined in 
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several different experiments. It was observed, that the accuracy of the 
interpolation results was sensitive to scaling. In general, the results with 
experimentally determined weighting of the normalized terms performed 
better than the version of the scheme that was neither normalized nor 
weighted. The available time for this investigation did not allow for ex¬ 
tensive testing of different weighting values. Probably there is room for 
improvement of the method's accuracy by optimizing the weighting. The 
mean averaged absolute difference between observed and interpolated 
sound speed values in the experiments was on the order of 0.5 m s. This 
seems like a relatively good preliminary result. However, to detect the 
expected greenhouse signal, it will be necessary to determine the variable 
acoustic travel times with at least a 10 ms accuracy (Munk and Forbes, 
1989). For the acoustic transmission simulations in this thesis the modi¬ 
fied multiquadric scheme would therefore not be sufficiently accurate. 
The modified MQ scheme performed better, when interpolants were 
sought at levels that were also used as grid point depth levels. Possibly 
this is due to the fact, that the additional condition that the sum of the 
MQ coefficients is equal to zero was not applied here. Effectuating this 
condition in the algorithm could improve the method's accuracy in future 


investigations. 




The method retrieves grid point data exactly within necessary preci¬ 
sion. This reflects the theory behind the method. The size of the area 
has definite effects on the degree to which the entries in the MQ coeffi¬ 
cient matrix are distinct. Kansa (1990) discussed the influence of this 
distinctness on the goodness of fit of the method. Future research should 
investigate this matter further. It is suggested to perform testing for larger 
areas and more vertical levels. Smaller values for d in general perform 
better than larger values. The experiments suggested that there exists an 
optimum value for <5; for this case equal to 10" 4 . Kansa's suggestion of 
varying the values of S with basis function number was not further ex¬ 
plored here. Future investigations should not overlook this suggestion. 
The goodness of fit of the method generally improved when more grid 
points are used. 

The setup and available time for these experiments did not allow for 
more general and definitive conclusions. The modified multiquadric 
method seems a potentially useful method for application to global ocean 
acoustics experiments, especially because corrections for the ellipticity of 
the earth can be incorporated rather easily. However, apart from this in¬ 
vestigation, the method has not yet been used in this field at all and more 
research is needed. The accuracy of the scheme is, as yet, not sufficient 
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for application to the global cross-basin range acoustic transmission 
simulations, ft was therefore decided to select the EOF method, which is 
discussed below. Nevertheless, improvements to the method along the 
lines suggested could “upgrade” the scheme for actual future use in 
acoustics projects. 

C. EMPIRICAL ORTHOGONAL FUNCTIONS METHOD 

In the upgraded version of HARPO, which is used in this thesis, the 
interfacing method of choice involves the Empirical Orthogonal Fun 
tions (EOF) representation. In contrast to the experimental nature of the 
modified multiquadric scheme, the EOF method has been routinely ex¬ 
ploited and tested by the physical oceanic community (Newhall, ct al., 
1989). To begin the description of the EOF scheme, recall how a rec¬ 
tangular window is used to move in a discrete fashion with the ray being 
traced. Data inside a sub-region with boundaries defined by the window 
are extracted. Inside this sub-region, sound speed is expressed as a sum 
of the mean sound speed c(z), a function of depth only, and the perturbed 
speed <5c. In formula form: 

c{6,(j),z) = c(z) + dc{d,(f),z ) (17) 
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The perturbation is then expressed as a linear combination of empirical 
orthogonal functions (EOF's): 


X 

Sc(d, <j>, z) = a t (6, (z) 

/=! 


08 ) 


where N is at most the number of vertical layers in the ocean model 
(i.e., 20) , a,{6,4 >) are the horizontal coefficients and J(z) are the vertical 
EOF's. This modal representation has two immediate advantages 
(Newhall, et al., 1989): 

• in representing a three-dimensional function as a product of two- 
dimensional and one-dimensional functions, one is able to use simple 
2-D and 1-D spline routines to supply smooth functions for the field 
at all points; 

• the EOF's are the most efficient energy representation available, so 
that one needs the fewest number of modes to account for a given 
percentage of observed variance. 


The latter advantage allows one to optimally truncate the mode sum, 
when the higher modes contain little energy. Preliminary investigations 
showed, however, that the cross-basin acoustic ray paths are very sensitive 
to vertical resolution. Ray trajectories change significantly as the number 
of modes that are retained is changed. Consequently, the maximum 
number of modes, 20, was used for the analysis presented in this thesis. 
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The next step involves obtaining the mean profile, the EOF's, their 
coefficients and their energy. A layered mean sound speed profile c[c k ) 
is calculated by averaging all the layered profiles at each horizontal 
gridpoint within the window, such that the mean at the k‘ h vertical layer 
is: 


?(--*> = ir 1 c a 0 j- o’) 

/= 1 

where M is the number of horizontal gridpoints. In the next step, em¬ 
pirical orthogonal functions of the perturbation profiles are constructed 
by solving the following eigenvalue / eigenvector problem: 

A r A CM 

Thus, the EOF's, f, are the eigenfunctions of the covariance matrix 
A r A and /, are their associated energies. To solve this equation singular 
value decomposition is used, i.e., set 

A r A=UAV 7 (21) 

Using this form, the columns of the U matrix are the f and A gives the 
mode energies. Next, the horizontal coefficients of the EOF's are 
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computed by projecting the layered sound speed perturbation profiles on 
the layered EOF's, i.e., 

N 

<*i 4>) = S z nVi ( z «) = 1 X..;N (22) 

n= 1 

where i is the mode index and n is the layer index. Cubic splines are then 
fitted to the mean profile as well as to the EOF's and bi-cubic splines are 
fitted to the EOF coefficients. This procedure will result in a smooth and 
continuous sound speed field, as required by the HARPO code. 

After the ray is traced to a point near the edge of the window, the 
window is moved to center at that position. New sub-region data are then 
extracted and the process is repeated. 
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IV. ACOUSTIC RAY VARIABILITY ANALYSIS 


In this chapter the procedure for and the results of the analysis of 
acoustic ray variability under the influence of meso- and gyre scale ocean 
temperature fluctuations are presented and discussed. Emphasis is placed 
on examining the temporal and spatial variability of the wave front travel 
times and azimuthal arrival angles of individual rays, as well as the lo¬ 
cations along the west coast of the United States being insonified by an 
envelope of reliable rays launched from a Heard Island sound source and 
the temporal variability of this ray envelope. The procedure for the 
analysis along with some results concerning the determination of reliable 
rays will be discussed first. 

A. PROCEDURE 

1. Programming Code Modifications 

For the ray variability analysis raytracing was performed using the 
upgraded three-dimensional Hamiltonian raytracing code HARPO. The 
code was successively interfaced with various output data sets from the 
Semtner-Chervin global ocean eddy resolving general circulation model, 
using a empirical orthogonal functions (EOF) method, as described in 
Chapter III, Section C. 
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Apart from the modifications involved in the interfacing proce¬ 
dure, two additional changes were be made in the HARPO code. First, 
the original HARPO output data set did not contain information about 
the ray point coordinates r, 0 and d or the sound speed and sound speed 
gradient along the ray trajectory. Since such information was necessary 
for this analysis Ehret (1990) modified the code to produce an extra out¬ 
put data set, i.e., a “rayset”, for each HARPO run. containing this sup¬ 
plementary raypath information. In this data set, output values were 
stored eveiy five integration steps. The rayset-files were subsequently used 
in all further calculations and processing. 

The second change corresponds to approximating the coastline by 
a straight line in the horizontal plane. This artificial coastline acts as a 
consistent terminator for the rays. The line was chosen so as to conform 
as close as possible to the geometry of the Semtner-Chervin ocean model 
in the vicinity of Monterey Bay, while maintaining its linear character, 
necessary to simplify the interpretation of the raytracing results. In 
radians this artificial coastline was defined as: 

0=1.25 0- 4.258603376 (23) 
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2. Selecting Reliable Rays 

With the modified HARPO code, rays were traced from the Heard 
Island sound source location using the 900-day mean field ocean model 
output data set. Subsequently rays were traced using the 11 consecutive 
instantaneous fields (hereafter to be referred to as the FI through Fll 
fields), extracted from the ocean model output data at 30-day simulation 
time intervals. Because cross-basin rays that encounter frequent bottom 
interactions suffer significant losses, the corresponding arrivals at listening 
stations close to the west coast of the United States are unreliable. It is 
therefore necessary to determine the limiting elevation angles between 
which the “reliable” rays are confined. A reliable ray is defined here as 
a ray that is not impeded by land masses and that has fewer than five 
bottom interactions along its trajectory. To determine the limiting launch 
elevation angles, rays were traced using a 133° launch azimuth angle in 
the mean field setup. The launch elevation angle was varied between 
+ 4° and -8°, with a 2° interval. It was found that the 0° and -2° launch 
elevation angles both produced reliable acoustic paths, while the rays with 
positive or more negative elevation angles invariably became frequently 
bottom interacting (see for example Figure 5). Since a 0° launch ele¬ 
vation angle requires a nearly infinitely small integration step size in the 
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raytracing code, which can be troublesome under some circumstances, it 
was decided to carry out all further raytracing using a 2° downward 
launch elevation angle, representative of all reliable rays. These reliable 
rays will probably constitute the latest arrivals, making up the trailing 
edge of the observed arrival pattern at the receiver. Rays outside the 
envelope, that interact frequently with the bottom, will generally traverse 
higher sound speed water and will therefore arrive earlier. 

Next, rays were traced with the representative -2° elevation angle 
and varying launch azimuth angles, to determine the limiting azimuthal 
angles, between which rays would not be impeded by land masses. For 
the 900-day mean field ocean model output data set the launch azimuth 
angles were varied with a half degree increment, while for the FI through 
Fll fields an increment of one degree was used. For all fields it was 
found that only rays within an horizontal envelope between launch 
azimuth angles of 133° and 136° made it all the way to the west coast of 
the United States. Rays with smaller launch azimuth angles invariably 
ran into the Campbell Plateau, south of New Zealand (see Figure 3 on 
page 19), while rays with launch azimuth angles larger than 136° were 
impeded by Antarctica. Recapitulating, it was shown that all reliable rays 
fell within a ray envelope, that was limited vertically by elevation angles 
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Figure 5. A typical ray path for rays launched with elevation angles smaller than 
-2°: Vertical section along the ray trajectory for a ray with launch 
azimuth of 133° and launch elevation of -4°, using 900-day mean field 
ocean model output data set for a range from 0 to 250 km. This illus¬ 
trates frequent bottom interactions that rays with elevation angles less 
than -2° encounter.. 

of 0° and -2° and horizontally by azimuth angles 133° and 136°. The 
horizontal trajectories for all reliable rays within this envelope for the 
mean field are displayed in Figure 6. 


39 






70 90 1 10 130 150 170 190 210 230 

E. LONGITUDE (DEG. I 


Figure 6. Horizontal trajectories of reliable rays within the ray envelope for the 
900-day mean ocean model data field: Rays with the following launch 
azimuth angles arc displayed, from left to right at the source location: 
133, 133.5, 134, 134.5, 134.9, 135.5 and 136 degrees. 

Although the reliable ray envelopes for all fields were initially the 
same in terms of limiting launch elevation and azimuth angles, the 
acoustic paths that define the envelope at ranges away from the source 
differ from field to field under the influence of the modeled meso- and 
gyre scale ocean temperature fluctuations. 
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In the ray variability analysis emphasis was placed on those reliable 
rays in the described envelope, that made it from the proposed Heard 
Island sound source location (53.2°S / 73.1°E) through the Indian Ocean, 
passing south of New Zealand and through the Pacific Ocean to the west 
coast of the United States. Additionally the possibility of reliable propa¬ 
gation paths from the Heard Island location through the Indian Ocean, 
the Tasman Sea, north of New Zealand and the Pacific Ocean (to Coos 
Bay, Oregon) was also examined. In this case, using again a 2° down¬ 
ward launch elevation angle in the mean field setup, the azimuthal launch 
angles were varied between 111° and 119°. It was found that the ray 
launched with an azimuthal angle of 111 0 ran into Australia, while the 
119° ray was impeded by the northwest coast of New Zealand. All of the 
rays with azimuthal launch angles in between these limits have excessive 
bottom interactions in the Tasman Sea, due to its shallow bathymetry and 
the downward refracting nature of local sound speed profiles. A vertical 
section along a typical Heard Island to Coos Bay ray trajectory is dis¬ 
played in Figure 7. No reliable rays could thus be identified for this case. 

In carrying out the raytracing simulations it was found that the 
HARPO code did not handle all rays correctly. In some instances the 
integration process stopped somewhere along the ray trajectory Solving 
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Figure 7. A typical ray path passing through the Tasman Sea, showing frequent 
bottom interaction: Vertical section along the ray trajectory for a ray 
with launch azimuth of 114° and launch elevation of -2°, using the 
900-day mean field ocean model output data set and for ranges from 0 
to 18000 km. 

this programming error is deferred to future investigations. Whenever a 
ray calculation stopped before being finished, the ray was traced again, 
using a slight (generally ± 0.1°) deviation from its original azimuthal 


launch angle. 







For each of the rays traced two plots were produced using NCAR 
graphics: the First illustrating the horizontal ray trajectory in a cylindrical 
equidistant projection at the ocean surface; the second depicting a vertical 
section along the raypath. These plots, in combination with the HARPO 
output data sets, were studied carefully to determine whether the rays that 
were traced could be used reliably for the ray variability analysis. All rays 
that had excessive bottom interactions along their trajectory or were 
otherwise impeded by land masses, were disregarded and not used any 
further. 

3. Ray Arrival Positions 

To find an optimal location for a listening site on the west coast 
of the United States and to determine whether the shallow Erben and 
Fieberling seamounts, which were omitted in the modeled bathymetry, 
could possibly impede the arriving rays, it is necessary to know between 
what positions reliable rays are expected to arrive. For each selected, 
reliable ray a stop position on the predefined artificial coastline was 
computed by direct interpolation between two positions; one on either 
side of the coastline. The distribution of these arrival positions along the 
coastline will give a clear indication of the best listening array locations. 
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4. Wave Front Travel Times 

Interpreting the measured arrival pattern of the acoustic signals 
at the listening site will be greatly facilitated by the results of an analysis 
of the spatial and temporal variations in acoustic travel times. Travel 
times to the ray arrival position were directly interpolated from the rayset 
values. For estimation of wave front travel times to a fixed (receiver) 
depth, a travel time correction is applied to these values, rather than at¬ 
tempting to find eigenrays, that connect source and receiver. 

Whenever a ray does not arrive at the coastline at a depth exactly 
equal to the receiver depth, the wavefront travel time to a fixed depth 
will be either shorter or longer than the interpolated ray travel time, de¬ 
pending on the geometry of the particular ray arrival. There are four 
possible cases: the ray can arrive going upward or downward, at a depth 
below or above the receiver depth. It is possible to determine the arrival 
geometry for each ray from the rayset files. From the diagram in 
Figure 8 it can be seen that the wave front travel time correction is equal 
to: 


At — — 


!h sin(a)| 
- 


where c is the local sound speed. 


( 24 ) 
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As will be discussed below, the distance between the upper- and 
lower turning points of the rays, i.e., the ray vertical extent, at ranges 
larger than about 8000 km turned out to be an important factor in the 
travel time variability analysis. The double amplitude at a specified range 
of 12000 km was used for the classification of rays. Two different types 
of rays were found. This 12000 km range was chosen because, as will be 
discussed below, it assures that the rays are well away from the transition 
zone between latitudes 60°S and 40°S, where the largest sound speed 
variations occur. 

5. Azimuthal Arrival Angles 

To determine an optimal orientation of a possible horizontal re¬ 
ceiver array at the listening site on the west coast of the United States, it 
is necessary to knew angles of the arriving rays and the variability of these 
angles. For each of the reliable rays the azimuthal arrival angle with 
respect to north was therefore calculated using standard trigonometric 
formulas. 

B. RESULTS AND DISCUSSION 
1. Ray Parameter Variability 

Spatial variability of reliable rays within the ray envelopes is at¬ 
tributed to geographical variations in the sound channel axis depth and 


46 






sound channel strength. Different rays encounter different transverse and 
along-path sound speed gradients. Figure 9 displays a contour map of 
sound speed at the sound channel axis for the 900-day mean field and in 
Figure 10 the depths of these sound speed minima are contoured 1 . 

Rapid fluctuations of the sound channel axis depth and speed oc¬ 
cur between latitudes 60°S and 40°S, where the sound speed structure 
changes f.om a basically upward refracting half-channel regime in the 
cold Antarctic w'ater towards a deep sound channel regime with an asso¬ 
ciated shallow surface duct in the more northerly Pacific waters. In this 
transition zone the sound channel deepens about S00 m in just 10° 
equat^rward, while the axial sound speed increases by more than 20 m s 
over the same distance. Since the Heard Island sound source, at 
53.2°S / 73.1°E, is located in this transition zone, slight variations i.i 
azimuthal launch angle will cause large changes in sound speed along the 
ray trajectory. Away from the transition zone variations are much more 
gradual. 

In comparing the ray trajectories for the different fields, it was 
found that a general distinction could be made between two types of rays 

I Chiu, el al. (1990) found that maps of these two quantities, derived from the Scmtncr- 
Chervin ocean model output data used in this analysis, were very similar to those given by Munk 
and 1 orbes (1989). This finding fortified the confidence with which this ocean model was used 
here. 
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Figure 9. Sound speed at the sound channel axis for the 900-day mean field ocean 
model data set: Contours from 1450 to 1495 m/s, contour interval 
5 m/s. 

that were affected in a different way by the modeled sound speed 
structure: 

• although both types of rays start out in the upward refracting half¬ 
channel regime near Heard Island, the first general type of ray 
(Type I) is trapped closer to the surface than rays of the second type 
(Type II). Although the deep sound channel starts to appear at a 
range of about 5000 km the Type I rays remain trapped close to the 
surface by the shallow surface duct. Not until this surface duct 
weakens between 7500 and 8000 km ranges will the Type I rays enter 
the deep sound channel. The distance between the surface duct and 
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Figure 10. The depth of the sound channel axis for the 900-day mean Tield ocean 
model data set: Contours from 200 to 1800 m, contour interval 
200 m. 

the deep sound channel axis is rather large (about 1200 m), causing 
the rays to fluctuate about the axis with a large amplitude after the 
transition. The wide vertical separation between upper- and lower 
turning points (generally between 1000 and 2000 m) continues until 
the ray reaches the west coast of the United States. Since a large part 
of the ray trajectory traverses high sound speed water on either side 
of the sound channel axis, the travel times of this type of ray will 
generally be shorter than those of the second type of ray. Also the 
Type I rays will encounter larger sound speed variability than 
Type II rays because the upper parts of the ray trajectories lie within 
the oceanic thermocline where sound speed variability is largest. 
Type II rays are confined closer to the sound channel axis, well away 
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from the thermocline and will meet less pronounced sound speed 
variability. 

• the second type of ray is also trapped near the surface in the half¬ 
channel regime, but less close than Type 1 rays. At ranges around 
5000 km these rays encounter a stronger along-path sound speed 
gradient than Type I rays, which forces them to leave the surface 
duct and enter the deep sound channel below it. Around this range 
the deep sound channel axis depth gradually increases to about 
1000 m over a distance of 600 km. Because of this gradual transition 
Type II rays will fluctuate about the sound channel axis with a much 
smaller amplitude than Type I rays. Similar to Type 1 rays, the 
distance between upper- and lower turning points (for Type II rays 
generally less than 500 m) remains almost constant for ranges larger 
than about 8000 km. Because this type of ray is confined much closer 
to the sound channel axis, they will traverse lower sound speed water 
and in general will have longer travel times than rays of Type I. As 
discussed above, the sound speed variability which these rays en¬ 
counter along their trajectories will be less pronounced. 


All rays fall into either one of these two classifications. 
Figure 11 displays typical vertical sections along ray trajectories of a 
Type I and a Type II ray. For completeness, the horizontal trajectories 
of these rays are illustrated in Figure 12. 

In Tables 1 through 5 a summary of the parameters associated with 
each ray inside the ray envelopes for the various fields is presented. In 
addition to these tables, variability in ray parameters can be illustrated 
by means of scatter plots. In such diagrams, the values of certain ray 
parameters are plotted against arrival latitude along the straight coastline. 
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Figure 11 
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Typical vertical sections along the ray trajectories of a Type I and a 
Type II ray: Type I ray, launch azimuth angle 135°, F9 field (top); 
Type II ray, launch azimuth angle 133°, F9 field (bottom). Note that 
the modeled bottom bathymetry is also displayed. 
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Figure 12. Horizontal trajectories of typical rays for the F9 instantaneous ocean 
model data field: Type I ray with launch azimuth angle 135° on the 
right hand side; Type II y with launch azimuth angle 133° on the left 
hand side. 

The scatter plots for wave front travel times and azimuthal arrival angles 
are displayed in Figures 13 through 16. 

Looking first at the spatial variability in travel times from ray to 
ray within the ray envelopes for the 11 instantaneous fields, it can be seen 
that the average travel time spread in any of the envelopes is about 62 s. 
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Table 1. MEAN FIELD RAY PARAMETERS 


L.aunch 

azimuth 

(deg) 

Arrival E. 
longitude 
(deg) 

Arrival N. 
latitude 
(deg) 

Arrival an¬ 
gle w.r.t. 
North (deg) 

Wave front 
travel time 
(sec) 

Ray vertical 
extent at 

12k km 
range (m) 

Type 
of ray 

MEAN HELD 

132.0 

Impeded by Campbell Plateau 

132.5 

Impeded by Campbell Plateau 

133.0 

237.79 

36.77 

38.6631 

12101.1846 

420.0 

11 

133.5 

238.67 

35.67 

44.9658 

12019.6348 

1950.0 

I 

134.0 

237.92 

36.59 

44.9701 

12082.7373 

440.0 

11 

134.5 

238.08 

36.41 

36.8631 

12079.8018 

450.0 

II 

134.9 

239.07 

35.16 

44.9658 

12028.3662 

2220.0 

1 

135.5 

238.77 

35.53 

44.9886 

12029.4424 

2060.0 

I 

136.0 

237.94 

36.57 

34.9922 

12081.0947 

470.0 

11 

136.5 

Frequently bottom interacting near Antarctica 

137.0 

Ran into Antarctica 


Table 2. F1-F2 FIELD RAY PARAMETERS 


Launch 

azimuth 

(deg) 

Arrival E. 
longitude 
(deg) 

Arrival N. 
latitude 
(deg) 

Arrival an¬ 
gle w.r.t. 
North (deg) 

Wave front 
travel time 
(sec) 

Ray vertical 
extent at 

12k km 
range (m) 

Type 
of ray 

FI FIELD 

132.0 

Impeded by Campbell Plateau 

133.0 

238.19 

36.26 

39.8179 

12078.7139 

550.0 

11 

134.0 

238.04 

36.46 

44.9943 

12080.9873 

600.0 

II 

135.0 

238.45 

35.94 

45.0085 

12060.8691 

500.0 

II 

136.0 

238.01 

36.49 

35.5488 

12072.3945 

1000.0 

11 

137.0 

Ran into Antarctica 

F2 FIELD 

132.0 

Impeded by Campbell Plateau 

133.0 

238.52 

35.86 

45.0000 

12039.7881 

2000.0 

I 

134.0 

238.21 

36.23 

33.7268 

12056.3271 

1800.0 

1 

135.0 

238.67 

35.67 

36.8631 

12039.8086 

1960.0 

1 

136.0 

238.63 

35.71 

36.2552 

12017.5322 

1880.0 

1 

137.0 

Ran into Antarctica 

































































































































Table 3. F3-F6 FIELD RAY PARAMETERS 


Launch 

azimuth 

(deg) 

Arrival E. 
longitude 
(deg) 

Arrival N. 
latitude 
(deg) 

Arrival an¬ 
gle w.r.t. 
North (deg) 

Wave front 
travel time 
(sec) 

Ray vertical 
extent at 

12k km 
range (m) 

Type 
of ray 

F3 FIELD 

132.0 

Impeded by Campbell Plateau 

133.0 

238.43 

35.96 

36.8631 

12069.2793 

300.0 

II 

134.1 

238.17 

36.30 

44.9943 

12049.6924 

1790.0 

1 

135.0 

238.61 

35.73 

36.8631 

12041.4619 

1880.0 

1 

136.0 

238.04 

36.45 

45.0341 

12053.3643 

1900.0 

1 

137.0 

Ran into Antarctica 

F4 FIELD 

132.0 

Impeded by Campbell Plateau 

133.1 

238.44 

35.96 

36.8631 

12078.0840 

500.0 

II 

134.0 

237.98 

36.53 

41.1853 

12084.0000 

290.0 

11 

135.0 

239.00 

35.25 

36.8631 

12021.1396 

1820.0 

1 

136.0 

237.58 

37.02 

33.6927 

12109.4238 

560.0 

II 

137.0 

Ran into Antarctica 

F5 FIELD 

132.0 

Impeded by Campbell Plateau 

133.5 

238.09 

36.38 

44.9943 

12097.1357 

590.0 

11 

134.0 

237.97 

36.53 

45.0043 

12089.1748 

600.0 

11 

135.0 

239.00 

35.25 

33.7268 

12039.5293 

1200.0 

1 

136.0 

237.69 

36.90 

36.8631 

12091.6689 

950.0 

II 

137.0 

Ran into Antarctica 

F6 FIELD 

132.0 

Impeded by Campbell Plateau 

133.0 

238.21 

36.24 

36.9068 

12074.8457 

570.0 

II 

134.0 

237.95 

36.57 

36.8631 

12090.8115 

600.0 

II 

135.0 

238.45 

35.93 

36.8631 

12115.6572 

1200.0 

II 

136.0 

l- 

237.92 

36.59 

38.6298 

12083.0791 

950.0 

11 


137.0 


Ran into Antarctica 


54 

































































































































Table 4. F7-F10 FIELD RAY PARAMETERS 


Launch 

azimuth 

(deg) 

Arrival E. 
longitude 
(deg) 

Arrival N. 
latitude 
(deg) 

Arrival an¬ 
gle w.r.t. 
North (deg) 

Wave front 
travel time 
(sec) 

Ray vertical 
extent at 

12k km 
range (m) 

Type 
of ray 

F7 FIELD 

132.0 

Impeded by Campbell Plateau 

133.5 

238.01 

36.48 

40.5949 

12085.8633 

600.0 

II 

134.0 

237.87 

36.66 

36.8631 

12096.8027 

400.0 

II 

135.0 

238.43 

35.96 

35.5229 

12067.9502 

290.0 

11 

136.0 

237.96 

36.55 

38.6631 

12103.9189 

740.0 

II 

137.0 

Ran into Antarctica 

F8 FIELD 

132.0 

Impeded by Campbell Plateau 

133.0 

238.16 

36.30 

36.8631 

12092.5117 

460.0 

II 

134.1 

238.44 

35.96 

41.1853 

12014.9453 

2000.0 

1 

135.0 

238.70 

35.62 

33.6927 

12051.8779 

830.0 

II 

136.0 

238.37 

36.04 

38.6665 

12079.1387 

510.0 

11 

137.0 

Ran into Antarctica 

F9 FIELD 

132.0 

Impeded by Campbell Plateau 

133.0 

238.19 

36.27 

44.9658 

12075.1787 

460.0 

II 

134.1 

237.81 

36.73 

45.0427 

12094.2061 

490.0 

II 

135.0 

239.55 

34.56 

44.9943 

11984.6299 

1810.0 

1 

136.1 

237.97 

36.54 

38.6665 

12093.0654 

610.0 

II 

137.0 

Ran into Antarctica 

FIO FIELD 

132.0 

Impeded by Campbell Plateau 

133.0 

238.91 

35.36 

44.9658 

12006.6113 

1900.0 

1 

134.0 

237.89 

36.65 

38.6648 

12055.5625 

1840.0 

1 

135.0 

239.15 

35.07 

36.8631 

12011.4775 

1820.0 

I 

136.0 

238.62 

35.73 

36.8631 

12000.8147 

2180.0 

I 

137.0 

Ran into Antarctica 


55 































































































































Table 5. Fll FIELD RAY PARAMETERS 


Launch 

azimuth 

(deg) 

Arrival E. 
longitude 
(deg) 

Arrival N. 
latitude 
(deg) 

Arrival an¬ 
gle w.r.t. 
North (deg) 

Wave front 
travel time 
(sec) 

Ray vertical 
extent at 

12k km 
range (m) 

Type 
of ray 

HI HELD | 

132.0 

Impeded by Campbell Plateau | 

133.0 

239.40 

34.74 

41.6392 

11964.6611 

2200.0 

I 

134.0 

238.55 

35.81 

39.2828 

12010.5859 

1920.0 

1 

135.1 

239.31 

34.86 

36.8849 

12005.1914 

1860.0 

I 

135.9 

238.16 

36.30 

44.9886 

12096.1309 

620.0 

II 

1 37.0 

Ran into Antarctica ] 


The maximum wave front travel time spread is 131.5 s in the Fll field 
envelope and a minimum spread of 20.1 s is found in the ray envelope 
of the FI field. The distinction between the Type 1 and Type 11 ray travel 
times is confirmed in Figure 14. The Type I rays generally have shorter 
travel times than the Type II rays. The travel times of the two groups have 
little overlap. 

Apart from this spatial variability, the wave front travel times also 
exhibit fluctuations that are due to variability of meso- and gyre scale 
ocean processes in time. A rough estimate of the temporal variability in 
travel time can be obtained by examining the travel time values found 
along lines of constant latitude. The scattered travel time values at a given 
latitude are an indication of the temporal rms fluctuations from field to 
field. A visual estimate of the rms temporal wave front travel time 












































N. LATITUDE (DEG. ) 


Figure 13. Wave front travel time scatter plot as a function of latitude, showing the 
distinction between different fields: Travel times labeled with M arc for 
mean field rays, labels 1-11 are for the Fl-Fll fields respectively. 











Figure 14. Wave front travel time scatter plot as a function of latitude, showing the 
distinction between Type I and Type II rajs: Travel times labeled with 
I are for Type 1 rays, label 11 applies to Type II rays. 











ARRIVAL ANGLE (DEG. 



Figure 15. Azimuthal arrival angle scatter plot as a function of latitude, showing the 
distinction between different fields: Arrival angles labeled with M arc 
for mean field rays, labels 1-11 arc for the FI-Fll fields respectively. 
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N. LATITUDE (DEG. ) 


Figure 16. Azimuthal arrival angle scatter plot as a function of latitude, showing the 
distinction between Type I and Type II rays: Arrival angles labeled 
with I arc for Type 1 rays, label II applies to Type II rays. 
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variability is on the order of ± 15 s. Comparing the travel times for the 
Type 1 and Type II rays in Figure 14, it is clear that the temporal rms 
fluctuations in this ray parameter for the Type II rays (about ± 10 s), is 
almost twice as small as the rms fluctuations for the Type I rays (about 
± 20 s). Apparently the temporal variability is larger for rays that trav¬ 
erse larger distances away from the sound channel axis, i.e., the Type I 
rays. This is due to the large temporal sound speed variability in the 
oceanic thermocline, which is traversed by the Type I rays only. The 
variability along the sound channel axis, where the Type II rays are con¬ 
fined is smaller. 

In a recent investigation Semtner and Chervin (1990) quantified 
the temporal meso- and gyre scale effects, as well as seasonal effects, on 
acoustic travel times along three stationary axial acoustic paths. They 
found that the trav„i time to San F rai'^o fluctuates about 0.5 s around 
a 360-dav running mean 2 under ti e ,u;!u:.' ce of mesoscale ocean vari¬ 
ability. One of their conclusions, base ! ..at .hese findings, states that the 
mesoscale effects on acoustic travel times are not large enough to obscure 
the anticipated signal of global change in the ocean. However, in their 
analysis they ignored path changes induced by the ocean variability. 

2 It should be noted that the interval at which data were extracted from the ocean model 
output in the Scmtner-Chcnin investigation was 3 simulation days, compared to the 30-day simu¬ 
lation time interval between the 1 ; 1 through I I1 fields used in this thesis 






Comparing their results with the results presented above, it is obvious that 
the temporal wave front travel time fluctuations found in this thesis are 
at least one order of magnitude larger than the ones found by Semtner 
and Chervin. Apparently the additional effects of path changes that were 
simulated in this thesis, but were not considered by Semtner and Chervin, 
are not negligible. Whether these additional effects due to changes in the 
acoustic paths are so random that they can be averaged out to allow for 
the detection of a greenhouse signal, cannot be answered within the scope 
of this thesis. Future investigations using more samples should address 
this problem. 

Next, the variability in azimuthal arrival angles is examined in an 
analogous manner. An indication of the spatial variability is obtained 
by looking at the spread between extreme values inside the different ray 
envelopes. On the average this spread is about 8.0°. A maximum occurs 
in the F5 field envelope with 11.3°, while the minimum of 1.8° is found 
inside the ray envelope of the F6 Field. The temporal rms variability is 
estimated from the scatter diagrams in Figure 15 and Figure 16. It can 
be seen that this rms variability has an approximate magnitude of 
+ 4.0°. These figures illustrate how most arrival angles are concentrated 
around values of 37° and 45°. Unlike the effects on travel times, no clear 
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cut distinction between the spatial and temporal effects on azimuthal ar¬ 
rival angles of Type I rays and on those of Type II rays could be made. 

Looking at the scatter plots that show the distinction between 
Type I and Type II rays, it is evident that Type I rays generally arrive 
more to the south on the coastline than the Type II rays. This is due to 
the fact that rays of Type I traverse a much larger part of the water col¬ 
umn than rays of Type II. Consequently, they will accumulate more 
southward horizontal refraction while crossing the circumpolar frontal 
system. 

2. Insonification Area and Envelope Variability 

Apart from the analysis of the parameters of individual rays, such 
as wave front travel times and azimuthal arrival angles, it is also necessary 
to define, compute and analyze quantities that describe the ray envelope 
which insonifies the California coast. The mean of the arrival positions 
of the rays within an envelope is an estimate of the location on which the 
sound energy is focused, i.e., the envelope focus. The standard deviation 
of these arrival positions is a relative measure of how much the rays inside 
the envelope diverge away from the position of this focal point. A smaller 
envelope divergence (i.e., the standard deviation) implies a stronger envel¬ 
ope focus. The envelope width is defined as the width of the area which 
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is insonified by a particular ray envelope. The envelope center is the center 
of the insonified area. In Table 6 below, an overview of these envelope 
parameters is presented. 


Table 6. RAY ENVELOPE PARAMETERS 


Field 

Envelope 
focus longi¬ 
tude (deg) 

Envelope 
focus N. 
latitude 
(deg) 

Envelope 
center lon¬ 
gitude 
(deg) 

Envelope 
center \. 
latitude 
(deg) 

Envelope 

divergence 

(km) 

Envelope 
width (km) 

F1 

238.17 

36.29 

238.23 

36.21 

29.0 

72.6 

F2 

238.51 

35.87 

238.44 

35.95 

29.5 

75.1 

13 

238.31 

36.11 

238.33 

36.09 

37.3 

95.3 

F4 

238.25 

36.19 

238.29 

36.13 

87.1 

233.9 

F5 

238.19 

36.27 

238.34 

36.08 

81.1 

217.2 

F6 

238.13 

36.33 

238.19 

36.26 

35.8 

87.8 

F7 

238.07 

36.42 

238.15 

36.31 

35.4 

92.2 

F8 

238.42 

35.98 

238.43 

35.96 

32.2 

90.0 

F9 

238.38 

36.02 

238.68 

35.65 

113.8 

287.1 

F10 

238.64 

35.70 

238.52 

35.86 

78.4 

208.5 

FI 1 

238.86 

35.43 

238.78 

35.52 

86.3 

206.6 

Average 

238.36 

36.06 

238.40 

36.00 

58.7 

151.5 

St. dev. 

0.23 

0.28 

0.19 

0.24 

29.3 

75.3 


Clearly the results in Table 6 illustrate which area along the 
California coast is insonified by the Heard Island sound source, as well 
as the temporal variability of the insonification envelope. Considering 
all 11 envelopes, it can be seen that the focus positions fluctuate around 
an average longitude of 238.36°E and an average latitude of 36.06°N, with 
a standard deviation in longitude of about 0.23° and in latitude of about 
0.28°. The centers of individual ray envelopes vary around a mean 
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longitude of 238.40°E and a mean latitude of 36.00°N, with slightly 
smaller standard deviations of 0.19° and 0.24° respectively. The focal 
points for the individual envelopes are invariably different from the 
envelope centers and can be located on either side of these centers. The 
average distance between focus and centers is 15.2 km, a maximum of 
50.1 km is found for the F9 envelope and a minimum of 2.4 km for the 
F8 envelope. The strength of the foci, as expressed by the envelope di¬ 
vergences, varies largely. The maximum divergence is 113.8 km, the min¬ 
imum is 29.0 km and the standard deviation over all envelopes is 
29.3 km. The width of the area along the California coast that is 
insonified by the individual envelopes varies from 72.6 km for the FI 
field to 287.1 km for the envelope in the F9 field. The mean width is 
151.5 km with a standard deviation of 75.3 km. Generally, a larger en¬ 
velope width will imply a larger envelope divergence. 

C. THE 1991 FEASIBILITY EXPERIMENT 

The last section in this chapter is devoted to translating the above 
findings to practical implications for the planned 1991 feasibility exper¬ 
iment. The simulations and subsequent analysis in this thesis demon¬ 
strated the existence of unimpeded reliable acoustic path cross-basin 
range ray trajectories from the Heard Island sound source location to the 
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west coast of the United States. It was found that reliable rays fall inside 
a 3° x 2° angular sector bounded by launch azimuth angles of 133° and 
136° and launch elevation angles of 0° and -2°. Rays launched outside 
this angular sector will be impeded by the Campbell Plateau, south of 
New Zealand, or Antarctica, or will interact with the bottom frequently, 
thus rendering them unreliable. 

The expected positions of ray arrivals in the vicinity of Monterey Bay 
were shown to fluctuate between 34.56°N j 239.55°E and 
37.02°N / 237.58°E, with a mean focal point at 36.06°N / 238.36°E. A 
permanent Monterey Bay listening station should optimally be located 
near this position. A ship that deploys a listening array during the 1991 
feasibility experiment should provide experimental verification of this 
finding. 

Figure 17 illustrates the area of insonification and acoustic path den¬ 
sity for all imulated rays in the FI through FI 1 fields. From this figure 
it can be seen that the arriving rays stay well clear of the shallow 
Fieberling seamount, in position 32.4°N / 232.2°E. The other seamount, 
Erben, in position 32.9°N / 221SE falls outside of the plot boundaries 
and is also located well away from the envelope of arriving rays. Both 
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seamounts are thus expected to have no impeding effect on the ray ar¬ 
rivals in the vicinity of Monterey Bay. 

The possibility of reliable ray paths to Coos Bay, Oregon, was shown 
to be very slim. Rays within a azimuthal launch angle envelope between 
111° and 119° invariably interacted excessively with the shallow 
bathymetry in the Tasman Sea. It will probably be very difficult to detect 
the corresponding lower energy arrivals at a listening site near the Oregon 
coast. 

The wave front travel time that can be expected is approximately 
12060 s, with a variation of about 20 to 131 s. Most probably the reliable 
rays will constitute the latest arrivals, making up the trailing edge of the 
observed incoming acoustic arrival pattern. Estimated temporal rms 
fluctuations of the trailing edge at the receiver are on the order of 
± 15 s. 

Finally, it must be noted that the analysis did not allow for geodesic 
veering of acoustic rays under the influence of the earth's ellipticity. The 
HARPO code worked in polar, spherical coordinates and no ellipticity 
corrections were applied. To find an initial estimate of the magnitude 
of the error introduced by this omission, great-circle as well as geodesic 
trajectories were calculated using programs obtained from the National 
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E. LONGITUDE (DEG.) 


re 17. Area of insonification and acoustic path density in the vicinity of 
Monterey Bay: All ray arrivals for the FI through Fll fields arc 
shown. Position of Ficbcrling seamount is marked with an asterisk; 
Frbcn seamount in position 32.9°N / 227.5°F, falls outside of the chart 
boundaries. 










Oceanic and Atmospheric Administration 3 . In Figure 18 portions of 
great-circle and geodesic paths in the vicinity of Monterey Bay are dis¬ 
played. These non-acoustic paths have initial azimuth angles of 133, 134, 
135 and 136 degrees at the proposed Heard Island source location. From 
this plot it can be seen, that the geodesic paths always arrive more to the 
north than the great-circle trajectories. The distance between arrival po¬ 
sitions near the coast is on the order of 22 km. Based on this finding, it 
is expected that acoustic ray arrivals in reality will deviate to the north 
from their arrival positions presented in this thesis by about that distance, 
due to geodesic veering. Also it can be seen that the non-acoustic 
geodesic and great-circle paths are impeded by the shallow Fieberling 
seamount. However, acoustic rays are affected by large sound speed gra¬ 
dients in the transition zone, resulting in a southward displacement of 
about 350 km. As was shown above, this displacement due to horizontal 
refraction is large enough for the acoustic rays to stay well clear of the 
seamount. 


3 It must be emphasized, that the trajectories computed ir. this manner do not represent 
acoustic ray paths, but are the mere results of applying standard geodetic formulas. 
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E. LONGITUDE (DEG.) 


Figure 18. Great-circle and geodesic paths near Monterey Bay: Initial azimuth 
angles at the Heard Island source location were 133, 134, 135 and 136 
degrees, from bottom to top. In each pair of paths the geodesic is lo¬ 
cated to the north of the great-circle path. 













V. CONCLUSIONS 


In this final chapter the results and conclusions of this study will be 
summarized. Some final remarks along with recommendations for future 
work will also be made. 

A. SUMMARY OF RESULTS 

1. Reliable Rays 

In this study reliable rays were defined as rays that have less than 
five bottom interactions along their trajectories and that are not impeded 
by land masses. It was shown that all reliable rays will be confined to a 
ray envelope, bounded by launch elevation angles of 0° and -2° and 
azimuthal launch angles of 133° and 136°. The envelope dimensions are 
invariant in time near the source, but can be quite variable at the 
California coastline. 

2. Area of Ir.scnification and Optimal Receiver Site Location 

One of the primary objectives of this study was to determine the 
locations along the west coast of the United States that will be insonified 
by the acoustic signals emitted by the Heard Island sound source. Using 
numerical models it was demonstrated that reliable, unimpeded acoustic 
paths over cross-basin ranges do exist. The expected arrival positions of 
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reliable rays launched at Heard Island fluctuate between 
34.56°N / 239.55°E and 37.02°N / 23 7.58°E, roughly from San Francisco 
to Santa Barbara on the California coast, with a mean focus at 
36.06°N / 238.36°E. This suggests that a receiver array should optimally 
be located at, or close to, this position. Standard deviations of the en¬ 
velope foci were about 0.25° and the centers of the insonification envel¬ 
opes in general do not coincide with these foci. The average strength of 
the focal points, as expressed by the envelope divergences, was 58.7 km, 
but with large variability. 

Ray arrivals near the Oregon coast are possible, but have a less 
reliable nature. Frequent bottom interactions in the relatively shallow 
Tasman Sea will deplete the acoustic energy of these rays so much that 
it will be hard, if not impossible, to detect the corresponding arrivals. 

3. Wave Front Travel Time Variability 

Estimated wave front travel times are on the order of 12060 s, al¬ 
though it was shown that spatial variations within an envelope can be as 
large as 131 s, depending on the envelope width. Estimated temporal rms 
fluctuations have a magnitude of about ± 15 s. 

An important finding is that a clear distinction can be made be¬ 
tween two types of rays with different characteristics. All rays that were 
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examined fell into either one of these categories. Type I rays, which are 
trapped close to the surface up to approximately 8000 km ranges and 
exhibit a large ray amplitude away from the transition zone, were shown 
to arrive earlier than Type II rays. Rays of the latter type enter the deep 
sound channel sooner, at ranges of about 5000 km and will have a much 
smaller ray amplitude away from the transition zone. Consequently, wave 
front travel times for this type of rays are longer than for the Type I rays. 
Temporal variability was found to be larger for the Type I rays than for 
the Type II rays. This is due to the fact that Type I rays traverse larger 
distances away from the sound channel axis, partly through the oceanic 
thermocline, where sound speed fluctuations are expected to be more 
pronounced. Type II rays are confined much closer to the axis, away 
from the thermocline. Because of lesser contamination by the meso- and 
gyre scale variability Type II rays appear to be better suited for the de¬ 
tection of climate change in the deep ocean. However, only Type I rays 
contain information about ocean temperature fluctuations in the upper 
part of the water column. 

4. Azimuthal Arrival Angles 

Estimated azimuthal arrival angles vary between approximately 
33.7° and 45.0°, with an average of about 39.0°. A permanent (horizontal) 




receiver array should optimally be oriented in a direction close to this 
average. Expected temporal rms fluctuations in arrival angles are on the 
order of ± 4.0°. An optimal receiver beam-width should be on the same 
order of magnitude. 

5. Shallow Seamounts 

It was found that the shallow Fieberling and Erben seamounts, 
close to the California coast, will not impede the acoustic transmissions. 
Both are located well to the north of the incoming acoustic rays. 

6. Interfacing Procedure 

To aid in the continuous effort to examine whether existing 
ocean'acoustic interfacing procedures and methods in ongoing acoustic 
computer simulation studies can be improved upon in terms of accuracy 
and efficiency, two alternative methods to interpolate gridded sound 
speed data were examined. The results indicate that the EOF method 
developed by Newhall, et al. (1986) is better suited for this kind of study. 

7. Naval Operations 

In the context of antisubmarine warfare operations an important 
finding is that reliable acoustic paths over cross-basin ranges can exist. 
Depending on source levels and frequencies, detection of submarines over 
such long ranges might be possible. Furthermore it was shown how 
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several important acoustic parameters are expected to vary under the 
influence of ocean it -so- and gyre scale ocean variability This knowledge 
about the variability in space and time of the acoustic wave fields that 
propagate through oceanic fronts and eddies may aid in future design and 
use of sonar and in the interpretation of the results in an operational 
setup. Finally, this thesis has demonstrated the feasibility of coupling 
global oceanic and acoustic models. Such integration can be applied to 
future on-board, long range acoustL propagation prediction systems in 
support of antisubmarine warfare operations. 

B. FINAL REMARKS AND RECOMMENDATIONS 

In this last section some Final comments will be maue. The scope and 
methods of the research imposed some ’limitations, which must be kept 
in mind when looking at the results. 

First, only a limited number of rays (generally aoout seven) were 
traced for each of the ocean model output data fields. The fact that the 
ruytracing had to be carried out in a small UNIX workstation environ¬ 
ment, in which each t ,tracing run took about 12-14 wall-clock hours, 
limited the number of runs that could be carried out in the available time. 
Results produced by applying statistical formulas, that inherently assume 
large samples, therefore only represent rough “estimates”. This thesis 
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should thus be viewed as an initial investigation, which must be expanded 
in the future. Carrying out a large number of additional raytracing runs, 
could greatly enhance the statistical significance of the results. Future 
availability of a supercomputer, on which these additional runs can be 
carried out simultaneously in a shorter time frame, is highly desirable for 
this kind of study. 

Whether the temporal variability in wave front travel times will or will 
not allow detection of the expected greenhouse signal in a permanent 
global climate monitoring system, depends on whether this variability can 
be reduced through averaging. Significantly more raytracing runs, to 
produce a long time series, could be helpful in solving this problem. 

Secondly, it was shown that rays are very sensitive to small variations 
in the sound speed structure. This is partly due to the fact that geometric 
raytracing is a high frequency approximation. To complement the results 
presented here and to aid in interpreting them, analyses using alternative 
techniques that take into account the finite frequency effects, such as 
dispersion and diffraction, will be necessary 4 . 

Thirdly, ocean currents were ignored in the setup of the analysis. The 
cumulative effects of currents are not necessarily insignificant, due to the 

4 Examples of such “full*wave** methods arc normal mode and parabolic equation analyses 
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very long path length under consideration. The same is true for seasonal 
effects. The ocean model output data sets used here, were obtained from 
model runs with seasonally constant forcing 5 . Future investigations 
should incorporate ocean currents, as well as seasonal forcing, to quantify 
the effects on the acoustic parameters discussed in this thesis. 

Another limiting factor results from the fact that the grid resolution 
of the Semtner-Chervin global ocean general circulation model is just 
marginally enough to resolve mesoscale eddies and fronts. When the next 
generation of supercomputers becomes available, the grid resolution 
should be increased to resolve smaller scale ocean fluctuations. Future 
acoustic simulations can then assess the impact of this smaller scale 
variability by using model output data fields with the enhanced grid 
resolution. 

Also, it was noted that this analysis did not allow for geodesic veering 
of acoustic rays under the influence of the earth's ellipticity. Future in¬ 
vestigations should incorporate ellipticity corrections in the raytracing 
code. 


5 Recently seasonal forcing was incorporated in the ocean model. Output over ten model 
years, extracted at 3 day simulation time interval, is presently available at the National Center for 
Atmospheric Research (Semtner, 1990) 
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